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Оптимизация замкнутых маршрутов 
на транспортной сети 


Предложен точный алгоритм ветвей и границ для решения замкнутой общей задачи коммивояжера. 
Среди решений с одинаковой стоимостью выбирается то, что содержит в себе наименьшее количество 
ребер. 


Постановка задачи 


В статье рассматривается задача нахождения замкнутого маршрута минимальной 
стоимости, который проходит по всем пунктам транспортной сети. 


Транспортная сеть представлена связным взвешенным графом Н =(,0 ), где 
Г — множество вершин, И -— множество ребер. Пункту Г сети соответствует верши- 
на 1еЕУТ графа Н, 


И =и, а отрезку дорожного полотна, связывающего два соседних 
пункта Ги /, отвечает ребро {1 Л ЕЦ . Каждому ребру {1 7 приписан вес 4, м 


равный стоимости проезда или расстоянию между пунктами Ги 7; № — множество 
действительных неотрицательных чисел. Предполагается, что 4, =Ч;. 


Требуется построить замкнутый маршрут, проходящий по всем вершинам графа 
Н иимеющий минимальную сумму весов ребер (минимальную стоимость). 

Сформулированная задача получила название общей задачи коммивояжера (ОЗК) [1], 
а ее решение — маршрутом коммивояжера. 

Маршрутом в графе Н называется чередующаяся последовательность вершин 
и ребер, а цепью — маршрут, все ребра которого различны. Простая цепь — это марш- 
рут, состоящий из различных вершин. Замкнутый маршрут называется циклом, а цикл, 
в котором все вершины различны, — простым. Простой цикл, проходящий по всем 
вершинам графа Н точно один раз, называется гамильтоновым циклом, или обходом. 
ОЗК, в которой требуется найти обход с минимальной суммой весов ребер, известна 
как гамильтонова задача коммивояжера (ГЗК) [1], [2]. Не всякий граф содержит га- 
мильтонов цикл, то есть не каждый граф гамильтонов. Следовательно, ГЗК в отличие 
от ОЗК не всегда имеет решение. 

Исследованию ОЗК посвящено ограниченное число работ, в которых в основном 
обсуждаются подходы к построению алгоритмов ее решения. Вместе с тем проблема 
нахождения оптимальных гамильтоновых циклов в полном графе всесторонне изучена. 
В данной статье развиваются результаты из [1], [2], ориентированные на сведение ОЗК 
к известной метрической задаче коммивояжера (ЗК). 
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Обоснование и описание алгоритма 


Известно, что ГЗК МР- трудна [3]. Этот же статус имеет ее частный случай, ко- 
торый состоит в определении гамильтоновости графа. Поэтому найти решение ГЗК 
или установить, что ее граф негамильтонов, можно только методами с экспоненциальной 
трудоемкостью вычислений. Так как множество циклов, проходящих по всем верши- 
нам связного графа, не пусто, то, очевидно, ОЗК поддается приближенному решению 
за полиномиальное время. Ниже изложен точный алгоритм решения ОЗК, построен- 
ный на ее сведении к метрической ЗК, которая состоит в нахождении гамильтонова 


цикла минимальной стоимости в полном графе Н, =(У,Е) с весами ребер {1, /} ‚ рав- 
ными кратчайшим расстояниям между вершинами Ги / графа Н = (У, ) 5] = рп. 

Каждому ребру {1, }еИ графа Н =(Т,И) поставим в соответствие множество 
А, простых цепей, в которых вершины Г и / являются концевыми. Множество А, 
содержит единственную цепь, состоящую из ребра {1 Л ‚ если в графе Н нет других 
простых цепей, соединяющих вершины Ги /. Выберем в А, цепь @,,, ребра которой 
в совокупности имеют наименьший вес Б(а,) . Назовем ее кратчайшей цепью между 
вершинами Ги /. В зависимости от значений 4, получим случай а), когда для всех 
ребер {1 7 ЕЦ выполняется неравенство 

4, <Б(а,), (1) 
и случай 6), когда оно нарушается хотя бы для одного ребра. Неравенство (1) харате- 
ризует весовые соотношения ребер в графе Н подобно неравенству треугольника 

а, <а,+аз, (2) 
которое должно выполняться для всех троек вершин 1, /, КЕЙ, 1*], /=К, полного 
взвешенного графа Н, =(У,Е) в метрической ЗК [1]. 

Найдем любым известным методом, например, алгоритмом Флойда — Уоршалла [1], 
кратчайшие цепи @, между всеми вершинами графа Н и определим их веса Б(а,) р 
Бл=Ьт [1 

Матрица [2(% )]. удовлетворяет неравенству (2) и вместе с матрицей [4 | за- 


дает полный граф Н„ =(У,Е), где каждое ребро {1,7} отвечает цепи @, в Н весом 


(а, (рис. 1. 


Рисунок 1 —а) Граф Н ; 6) полный граф Н„ 
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Пусть граф Н =(У,Ц ) гамильтонов, т, Т — оптимальные решения соответст- 
венно ГЗК и ОЗК, построенные в Н, О(т) и С(Т) - стоимости построенных реше- 
ний. Найдем обход о’ минимальной стоимости в полном графе Н т 

Утверждение 1. Если для всех ребер {1, /} е 0 гамильтонова графа Н = (7,0) вы- 
полняется неравенство (1), то Т=х, Ст) = Б(т) : 

Доказательство. Решение ГЗК т содержит п ребер множества ИЦ. В случае а) 


ребро {1,/} полного графа Н„ имеет вес Б(а,)=а если {1 ЛЕО, и Б(а,)>а; 


Й > 
иначе. Отсюда следует, что в Н„ точное решение ЗК о’ совпадает с т ‚а его стоимость 


не превышает стоимость произвольного маршрута, который содержит п и больше ре- 
бер, включая Т. 


Утверждение 2. Если хотя бы для одного ребра 57 гамильтонова графа 
Н = (У) неравенство (1) не выполняется, то ст) ы (с) : 

Доказательство. В случае 6) в полном графе Н.„ найдется по крайней мере од- 
но ребро {1, /}, которое в Н имеет вес 4, > Б(а, ). Если гамильтонов цикл о’ графа 
Но, проходит по ребру т Л ‚ то его стоимость больше стоимости соответствующего 


замкнутого маршрута в Н ‚, содержащего вместо ребра {1 Л =И цепь @,. Если цикл 


о’ не включает ребер, которые нарушают неравенство (1), то он совпадает с решени- 
ем ГЗК. 


Допустим, что граф Н = (Уи ) не гамильтонов. Тогда в любом случае а) или 6) 
стоимость оптимального решения ЗК о’ для полного графа Н,„ = (У, Е ) равна стоимос- 


ти решения ОЗК Т для графа Н . Поэтому маршрут Т можно найти в результате по- 
строения в графе Н., гамильтонова цикла о’и замены в нем каждого ребра {1 м } ЕЕ 


на цепь @,, из ребер множества 0. 


Приведенные рассуждения позволяют описать алгоритм решения ОЗК. 
50. Н = (Уи ) — связный взвешенный граф с множеством вершин Г, 


И =и,и 
множеством ребер И , [4, — матрица весов ребер графа Н ‚где 4, е К ‚если {1, }ЕЦ(, 
и а; = со иначе, 1, } = 1 п ь у — множество действительных неотрицательных чисел. 

51. Построить матрицу [2% | кратчайших цепей между всеми парами вершин 
графа Н и матрицу [2%], в которой элемент (1) равен весу Б(а,) цепи @,; 
матрицы @ ]. и [р (а, )]. определяют полный взвешенный граф Н, =(7, Е), где каж- 
дое ребро {1, /} заменяет цепь @,, в графе Н. 


52. В графе Н.„ найти обход минимальной стоимости о’ любым известным ме- 


тодом решения метрической ЗК. 
53. Построить оптимальное решение ОЗК Т в результате замены каждого реб- 


ра {7, /} гамильтонова цикла о’ на цепь @, графа Н. 
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Минимизация длины маршрута коммивояжера 


Количество ребер в маршруте называется его длиной [4]. В общем случае связ- 
ный граф Н содержит несколько замкнутых маршрутов коммивояжера минимальной 
стоимости, включающих разное число ребер. Представленный алгоритм строит какой- 
либо из них. Потребуем, чтобы длина маршрута коммивояжера была наименьшей. 

Для выполнения этого требования достаточно применить на шаге $2 классичес- 
кий алгоритм решения ЗК (метод Литтла) с дополнениями к правилу выбора вершины 
ветвления. 

В методе Литтла ветвление выполняется после приведения исходной матрицы 
стоимостей и представляет собой пошаговый процесс разбиения множества решений 
ЗК на непересекающиеся множества, изображаемые в бинарном дереве перебора в ви- 
де вершин [5]. Корню дерева решений соответствует вершина ©, которая обознача- 
ет множество всех допустимых решений ЗК. Она является начальной вершиной вет- 
вления. 


Пусть (а, — значение элемента (1, 7 ) приведенной матрицы стоимостей, 1, } = 


=1, п. Каждый элемент матрицы представлен дугой в ориентированном мультиграфе 
С, › где любые две вершины Г и / соединены парой дуг (1, /) и (1,1). 

Ветвление начинается с выбора дуги (,1 ) мультиграфа С’, которая, исходя из 
значений элементов приведенной матрицы, оценивается как наиболее подходящая для 
включения в искомое решение. Дуга (^,1 ) разбивает множество всех решений © на 
два подмножества. Решения одного подмножества включают дугу (к,1 ), а решения 
другого не содержат ее. Корень дерева © порождает две текущие активные вершины. 


Обозначим их соответственно ((%,/)°) и ((.1)-) . Для подмножеств ((%,/)°) и ((.1)-) 


вычисляются нижние границы стоимостей решений ЗК Ф((к,1 ) °) и 6((.1) -) ‚ в каж- 


дую из которых входит своя константа приведения [6]. Вершина с наименьшей ниж- 
ней границей выбирается в качестве вершины ветвления. 
По матрице, соответствующей выбранной вершине, определяется следующая ду- 


га (+, 5) ‚ инициирующая ветвление. Если Ф((к, 1) °) ыы ((*.1) -) ‚ то множество ((%,1) °) 
разбивается на подмножества ((%,) ‚(",5) °) и (4.1) (7.5) -) ‚ иначе множество (=. -) 


разбивается на подмножества (=. ‚(",5) -) и (в. (5) -). Вершина ветвления по- 


рождает пару текущих активных вершин, затем формируются соответствующие им 
приведенные матрицы, по которым находятся нижние границы для каждого из двух 
полученных подмножеств разбиения. Очередная вершина ветвления определяется сре- 
ди всех текущих активных вершин, представленных висячими вершинами в дереве 
перебора. Путь из корня в висячую вершину содержит часть допустимого решения ЗК, 
состоящую из включенных в процессе ветвления дуг. Алгоритм Литтла завершает 
работу, когда следующая вершина ветвления включает все дуги обхода в графе С. 
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Если на каком-либо шаге ветвления несколько активных текущих вершин полу- 
чают наименьшее значение для своих нижних границ, то метод Литтла не гарантирует 
решение ОЗК, минимальной длины. Поэтому, чтобы найти такое решение, в процес- 
се построения дерева перебора для каждой его вершины вместе с нижней границей 
стоимости определяется параметр, оценивающий длину частичного решения. 

В строящемся дереве перебора путь из корня в висячую вершину содержит ду- 


ги, включенные и не включенные в частичное решение. Обозначим Р; = {> а } 


множество дуг, включенных в частичное решение в результате построения пути из 


корня дерева в висячую вершину /. Каждая дуга (, У множества Р, в мультигра- 


фах С„ является дугой (*,>%, }, которой в графе Н соответствует кратчайшая цепь 


а ‚ соединяющая вершины р и 4 и содержащая Ц р,а) ребер. Тогда частичное ре- 


Ра › 
шение, полученное для вершины /, включает 


= 1». Дам (3) 
№ "М |=) 
ребер графа Н . В методе Литтла для каждой вершины /, порожденной на шаге вет- 


вления, к нижней границе стоимости ф р добавим оценку 
Г,=Г-Р,. (4) 


Из двух и более активных текущих вершин с наименьшим значением нижних гра- 
ниц, в качестве вершины ветвления выбирается та, для которой оценка (4) минимальна. 


Пусть в связном графе Н = (Уи ) Т - точное решение ОЗК с наименьшим чис- 
лом ребер |Т|. Тогда если |Т|= и, то из утверждений 1 и 2 следует, что граф Н - га- 
мильтонов. Если т > п, то граф Н либо не гамильтонов, либо такой гамильтонов граф, 


в котором, по крайней мере, для одного ребра неравенство (1) не выполняется. Дру- 
гими словами, алгоритм решения ОЗК Т минимальной длины является алгоритмом 
поиска решения ГЗК для графа Н в случае а). Для графа Н в случае 6) поиск реше- 
ния ГЗК выполняется алгоритмом, предложенным в [7]. 


Пример 


Найдем маршрут коммивояжера, содержащий минимальное число ребер для гра- 
фа, изображенного на рис. 1 а). После выполнения шагов 30 и 51 алгоритма точного 
решения СЗК получим матрицы: 


12345 12345 

о о 4 о 1 134 6 1 

а |= 2|® ®о® 3 2 а о 
р. [2(@,) 347 ® 10 5, 
4] © 3 хо ох 6 416 3 10 хх 5 

51 25 б® 5125 5 ® 
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1 и 3 4 5 

1 0 (1:52) (1,3) (1,5,2,4) (1,5) 

[а | = 20" РЭО 0 (2, 5;3) (2, 4) (2, 5) 
У [5 3 
4 
5 


(3.0 (3, 5, 2) 0 (3:5:2.4):. (353) 
(4.2.5.1): 42) 3255 0 (4,2, 5) 
(5, 1) (5, 2) (5, 3) (5, 2, 4) 0 


ша -ь-Ф 
шышоышньъь 
эр = = = 


На шаге $2 построим обход минимальной стоимости для полного графа Н „, пред- 
ставленного на рис. 1 6). Построение выполним модификацией метода Литтла, мини- 
мизирующей число ребер в решении ОЗК Т. На рис. 2 представлено дерево перебора, 
по которому находится обход минимальной стоимости в графе Н„ и маршрут ком- 
мивояжера минимальной длины в графе Н. 

15, 0 


4) © вы 18,0 
17,0 а 
(1,3) (1,3) (4,2 (4,2 т 
18,0 18,0 8: 160) 9) 
(.5 (.5 в [5 
О о (1,5) и (12) п, 
от 69 [606 5 / 69 09/9, 
20,2 со 20, 22, : 506 
9 © (9 09 а" д (9 
(4,3) | (4,3) 2.35) (5)° (1,4) (1.4) 
20,2 
вос. О] 65 в © 
(3,5) | 20, 0 22,0 
и (3,5) | [(,5) ы 
в) 6 (9 9 [о 63| |@3) 
56 ‚© ® 20, 15) (657 
(4,5) | |(4,5) (4,5) | | (4,5) 20,2 (5,4) 


20,2 20 20, 87 
20,2 20,2 (41) (40) 
(89) (39 


Рисунок 2 — Дерево поиска модифицированного алгоритма Литтла 
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Каждой вершине дерева приписано значение нижней границы стоимости частич- 
ного решения и число, вычисляемое по формулам (3) и (4). Граф Н содержит несколь- 
ко оптимальных маршрутов коммивояжера равной длины. Например, ОЗК находится 
из пути от корня © дерева до вершины 37. Этот путь включает дуги (4, 2), (1,5), (3, 1), 
(2,3), (5,4), образуя в полном графе Н„ цикл (1, 5, 4, 2, 3, 1) стоимостью 20. 

На шаге 53 выполним замену ребер {5, 4} и {2, 3} соответственно на цепи 
ад = (5,2, 4) и а. = (2, 5, 3). В результате получим решение ОЗК Т =(1, 5, 2, 4, 2, 5, 
3, 1), состоящее из 7 ребер. Другого решения стоимостью 20 и с меньшим числом ре- 
бер в графе Н нет. 


Заключение 


В данной статье нашли развитие изложенные в [1], [2] идеи, которые лежат в ос- 
нове метода построения замкнутого маршрута минимальной стоимости, проходящего 
по всем пунктам транспортной сети. Задача построения такого маршрута называется 
общей задачей коммивояжера (ОЗК). Моделью транспортной сети является связный 
неполный граф с произвольно заданными неотрицательными весами (стоимостями) 
ребер. Показано, что если в графе содержится гамильтонов цикл, т.е. простой цикл, вклю- 
чающий каждую вершину точно один раз, то его стоимость не меньше стоимости реше- 
ния ОЗК. 

Предложен точный алгоритм решения ОЗК, выполняемый в две стадии. На пер- 
вой стадии граф транспортной сети преобразуется в полный граф с весами ребер, огра- 


ниченными неравенствами треугольника. Каждому ребру т, Л полного графа ставится 


в соответствие цепь минимальной стоимости, соединяющая вершины Ги /. На второй 


стадии применяется классический алгоритм ветвей и границ для решения метрической 
симметричной задачи коммивояжера, дополненный способом минимизации числа ребер 
в искомом маршруте. 
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А.В. Панишев, А.Ю. Левченко, О.Б. Маций 

ОптимзащЯя замкнених маршрут!в на транспортнй мереж! 

Запропоновано точний алгоритм плок та меж для розв’язку замкнено] загально! задач! ком1вояжера. 
Серед розв’язк1в з однаковою варт1стю обираеться той, що мстить найменшу юльксть ребер. 


А.Г. Раш5йеу, А. Уи. Геусйенко, О.В. Мау 

Тгапзрогё Мебуогк?”$ С1озеа У/’аЩ$ Орйпи7айоп 

Те агие оНег$ ехас{ БгапсВ ап Боипа$ ао ил оф с1о5е4 соттоп Соттегсла| Тгауег ТазК зоайоп. 
ТЬы$ а[еоги фил з6]ес{5 Фе зошоп, \УМсВ сошатз пишита| пипбег оЁ едез Бебмееп еда1 с0${ зо[аНопз. 
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